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As photons do not interact with each other, it is interesting to ask whether photonic systems 
can be modified to exhibit the phases characteristic of strongly coupled many-body systems. We 
demonstrate how a Mott insulator type of phase of excitations can arise in an array of coupled elec- 
tromagnetic cavities, each of which is coupled resonantly to a single two level system (atom/quantum 
dot/Cooper pair) and can be individually addressed from outside. In the Mott phase each atom- 
cavity system has the same integral number of net polaritonic (atomic plus photonic) excitations 
with photon blockade providing the required repulsion between the excitations in each site. Detun- 
ing the atomic and photonic frequencies suppresses this effect and induces a transition to a photonic 
superfluid. We also show that for zero detuning, the system can simulate the dynamics of many 
body spin systems. 



Introduction: Strongly coupled many-body systems de- 
scribed by the Bose-Hubbard models [l[ exhibit Mott 
insulating phases whose realization in optical lattices 
0) S 0] have opened varied possibilities for the simula- 
tion of many body physics [5] . Can there be another en- 
gineered quantum many-body system which displays such 
phases? This will be especially interesting if the strengths 
of this system are "complementary" to that of optical lat- 
tices - for example if it allowed the co- existence of acces- 
sibility to individual constituents of a many-body system 
and a strong interaction between them, or if it allowed the 
simulation of arbitrary networks rather than those deriv- 
able from superposing lattices. Particularly arresting will 
be to find such phases by minimally modifying a sys- 
tem of photons which, by being non-interacting, are un- 
likely candidates for the studies of many-body phenom- 
ena. Here we propose such a system consisting of coupled 
electromagnetic cavitie s d, 3, [ 9|. [lol 111, doped with 
single two level systems [12 
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Using the nonlinearity generated from the corresponding 



photon blockade effect [21 



22l | , we show the possibility of 
observing an insulator phase of total (atomic+photonic 
or simply polaritonic) excitations and its transition to 
a superfluid of photons. Compared to the optical lat- 
tices case, the different Hubbard-like model which de- 
scribes the system involves neither purely bosonic nor 
purely fermionic entities, the transition from insulator 
to superfluid is also accompanied by a transition of the 
excitations from polaritonic to photonic, and that exci- 
tations, rather than physical particles such as atoms are 
involved. In addition, the possibility to simulating the 
dynamics of an XY spin chain-with individual spin ma- 
nipulation, using a mechanism different from that used 
in optical lattices [5] is suggested. 

System Description: Assume a chain of N coupled 
cavities. A realization of this has been studied in 
structures known as a coupled resonator optical waveg- 



uide (CROWs) or couple cavities waveguides(CCW) in 
photonic crystals d, M, flol- [23t . in tabered fiber cou- 
pled toroidal microcavies ll4l and coupled superconduct- 

This has already 



ing microwave resonators |12l 1171 . [1 
stimulated proposals to use such systems to implement 
optical and cluster state quantum computingfi^, [iH l. 
produce entangled photons ,_26] and study interacting 
polaritons 27| . We will describe the system dynamics 
using the operators corresponding to the localized eigen- 
modes (Wannier functions), aj,(afc). The Hamiltonian is 

given hyH = Y.k=i '^da\ak + Z^fcLi A{a\ak+i + H.C.) 
and corresponds to a series quantum harmonic oscilla- 
tors coupled through hopping photons. The photon fre- 
quency and hopping rate is u)d and A respectively and 
no nonlinearity is present yet. Assume now that the 
cavities are doped with two level systems (atoms/ quan- 
tum dots/superconducting qubits) and \g)k and \e)k their 
ground and excited states at site k. The Hamiltonian de- 
scribing the system is the sum of three terms. H^^'^'^ the 
Hamiltonian for the free light and dopant parts, iJ™* the 
Hamiltonian describing the internal coupling of the pho- 
ton and dopant in a specific cavity and H^°p for for the 
light hopping between cavities. 

N 

ff/^- =UJdJ2alak+iOoJ2\e)k{e\k (1) 

k=l k 
N 

W^' =gJ2{al\g),{eU + H.C.) (2) 

N 

jjhop =A^(4a,+i+ff.C) (3) 



k=l 



where g is the light atom coupling strength. The H^^'^^ -\- 
£ji7it pg^j-^ Qf f^i^Q Hamiltonian can be diagonalized in a ba- 
sis of mixed photonic and atomic excitations, called po- 
laritons. These polaritons, also known as dressed states. 
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FIG. 1: A series of coupled cavities coupled through light and 
the polaritonic energy levels for two neighbouring cavities. 

involve an equal mixture of photonic and atomic exci- 
tations and are defined on resonance by creation oper- 
ators P^^'"^^ = (|5,n)fc(g,0U ± |e,n - l)k{gMk)lV2, 
where |n)fc, \n — l)k and |0)fc denote n,n — 1 and pho- 
ton Fock states in the fcth cavity. The polaritons of 
the kth atom-cavity system are denoted as \n±)^. and 
given by \n±)k = {\9,n)k ±\e,n- l)k)/V2 with ener- 
gies = nojd ± g\/n and are also eigenstates of the 
the sum of the photonic and atomic excitations operator 
J^k = a\ak + |e)(e|fe with eigenvalue n{¥ig. 1). 

Polaritonic Mott State: We will now justify that the 
lowest energy states of the system consistent with a given 
number of net excitations per site (or filling factor) be- 
comes a Mott state of the net (polaritonic) excitations for 
integer values of the filling factor. To understand this, 
we rewrite the Hamiltonian (for A 0) in terms of the 
polaritonic operators (defined in the caption of Fig. 1) 
as 

N oo 
oo 

^n(c., + .9)p(-^-)tp(+.n)^ 

oo 

oo 

^5(^AI~.)P^"M^-)] + 

77 = 1 

N 

Aj2{aW+i+H.C). (4) 

fc=i 

The above implies (assuming the regime A << g\fn << 
cijfi) that the lowest energy state for a given number, say 
77, of net excitations at the fcth site would be the state 
\r\—)k (this is because |r?-|-) has a higher energy, but same 
net excitation 77). Thus one need only consider the first, 
third and last lines of the above Hamiltonian Yl for de- 



termining the lowest energy states. The first line cor- 
responds to a linear spectrum, equivalent to that of a 
harmonic oscillator of frequency — 9- If only that 
part was present in the Hamiltonian, then it would not 
cost any extra energy to add an excitation (of frequency 
u!d — g) to a site already filled with one or more excita- 
tions, as opposed to an empty site. However, the term 
g{n — y/v^P'ir '"'^^ pIt '"'^ raises energies of uneven excita- 
tion distribution such as l(ri-l-l)— )fc|(n— 1)— ); among any 
two sites fc and I relative to the uniform excitation distri- 
bution \'n—)k\'n—)i among these sites. Thus the third line 
of the above Hamiltonian can be regarded as an effective, 
nonlinear "on-site" photonic repulsion, and leads to a 
Mott state of the net excitations per site being the ground 
state for commensurate filling. Reducing the strength of 
the effective nonlinearity, i.e., the blockade effect, should 
drive the system to the superfluid regime. For this, one 
should move the system away from the strong resonant 
interaction to a weaker dispersive regime. This could be 
done by Stark shifting and detuning (globally again) the 
atomic transitions from the cavity by an external field. 
The new detuned polaritons are not as well separated as 
before and their energies are merely shifts of the bare 
atomic and photonic ones by ±g^(n -I- I)/ 5 for the |e,n) 
and |(7, n ~\- \) respectively. In this case it costs no extra 
energy to add excitations (excite transitions to higher 
polaritons) in a single site, and the system moves to the 
superfluid regime. 

To justify the transition of the system from a Mott 
phase to a superfluid phase as detuning A = ujQ—ujd'is in- 
creased, we have performed a numerical simulation using 
3-7 sites using the Hamiltonian of Eq.(l)-(3)(numerical 
diagonalization of the complete Hamiltonian without any 
approximations ) jsj]. In the Mott phase the particle 
number per site is fixed and its variance is zero (every site 
is in a Fock state) . In such a phase, the expectation value 
of the destruction operator for the relevant particles-the 
order parameter is zero. In the traditional mean field 
(and thus necessarily approximate) picture, this expec- 
tation value becomes finite on transition to a superfluid, 
as a coherent superposition of different particle numbers 
is allowed to exist per site. However, our entire system is 
a "closed" system and there is no particle exchange with 
outside. Superfluid states are characterized by a fixed 
"total" number of particles in the three site system and 
the expectation of a destruction operator in any given 
site is zero even in the superfluid phase. Thus this ex- 
pectation value cannot be used as an order parameter for 
a quantum phase transition. Instead we use the variance 
of total number of excitations per site, the operator A/fc, 
in a given site (we choose the middle cavity, but any of 
the other cavities would do) to characterize the Mott to 
superfluid phase transition. This variance var{Afk) has 
been plotted in Figl2] as a function of log^o A for a filling 
factor of one net excitation per site. For this plot, we 
have taken the parameter ratio g/A = 10^ {g/^ = 10^ 




FIG. 2: The order parameter as a function of the detuning be- 
tween the hopping photon and the doped two level system(in 
logarithimic units of the matter-light coupling g). Simula- 
tions include results for 3-7 sites, with and without dissipa- 
tion due to spontaneous emission and cavity leakage. Close 
to resonance (0 < A < 10~^), where the photon blockade 
induced nonlinearity is maximum (and much larger than the 
hopping rate), the system is forced into a polaritonic Fock 
state with same integral number of excitations per site(order 
parameter zero-Mott insulator state). Detuning the system 
by applying external fields and inducing Stark shifts {A > g), 
weakens the blockade and leads to the appearance of differ- 
ent coherent superpositions of excitations per site( a photonic 
superfluid).The increase in number of sites leads to a sharper 
transition as expected. 



gives very similar results), with A varying from ~ W~^g 
to ^ g and cud,uio ~ lO^g. We have plotted both ideal 
graphs (if neither the atoms nor the cavity fields under- 
went any decay or decoherence) and also performed sim- 
ulations explicitly using decay of the atomic states and 
photonic states in the range of g/max{K, 7) ^ 10"^, where 
K and 7 are cavity and atomic decay rates. 

These decay rates are expected soon to be feasible in 
toroidal microcavity systems with atoms IJ] and arrays 
of coupled stripline microwave resonators, each interact- 
ing with a superconducting qubit [3]. For these simu- 
lations we have assumed that the experiment (of going 
from the Mott state to the superfluid state and back) 
takes place in a time-scale of 1/A so that the evolution 
of one ground state to the other and back is adiabatic. 
The simulations of the state with decay have been done 
using quantum jumps, and it is seen that there is still a 
large difference of var{Afk) between the Mott and super- 
fluid phases despite the decays. As expected the effect of 
dissipation reduces the final value of order parameter in 
the superfluid regime(population has been lost through 
decay) whereas in the Mott regime leads to the intro- 
duction of fluctations again due to population loss from 
the state. The Mott {var{Mk) = 0) to superfluid 
(var(J\fk) = 0.75) transition takes place over a finite vari- 
ation of A (because of the finiteness of our lattice) around 
lO^g and as expected becomes sharper as the number of 
sites is increased. 

In an experiment one would start in the resonant 



FIG. 3: Up: The probability of exciting the polaritons corre- 
sponding to one (thick solid line) or two (thick dashed line) 
excitations for the middle cavity which taken initially to be 
empty( in the resonant regime). In this case, as we started 
with fewer polaritons than number of cavities, oscillations 
occur for the single excitation polariton. The double one 
is never occupied due to the blockade efi'ect-Mott insulator 
phase(thick dashed line) . For the detuned case however, hop- 
ping of more than one polaritons is allowed (superfluid regime) 
which is evident as the higher polaritonic manifolds are being 
populated now(dashed line and thin solid line). 



(Mott) regime with all atom-cavity systems initially in 
their absolute ground (|g, 0)®*^) states and prepare the 
atom-cavity systems in the joint state |1— by apply- 
ing a global external laser tuned to this transition. This 
is the Mott state with the total (atomic -l-photonic) exci- 
tations operator A/fc having the value unity at each site. 
One would then Stark shift and detune (globally again) 
the atomic transitions from the cavity by an external 
field and observe the probability finding |1— ) and the 
predicted decrease of this probability (equivalent to the 
increase in our order parameter, the variance of A/fe) as 
the detuning is increased. For inferring the fluctuations 
in A4 of our system, it is basically suffices to check the 
population of the |1— ) state, as this is an eigenstate of 
A/fe. For this, a laser is applied which is of the right fre- 
quency to accomplish a cycling transition between |1— ) 
and another probe third level. Through monitoring its 
fluorescence, accurate state measurements can be made 
in the standard atomic state measurement way [2^. 

To strengthen our case, we have also calculated the 
probabilities of populating the lowest polaritonic state of 
the first and second manifolds |1— )fe, |2— )fe of a middle 
cavity out of an array of 3,5 and 7 cavities. The initial 
state is the polaritonic state \l—)k excited at the right 
and the left cavities, with the middle cavity being in the 
ground state \g,0)k- Figure [3] shows as expected from 
our discussion above, on resonance the photon blockade 
is preventing any excitation to any state higher than the 
first(Mott insulator phase). However, by simply vary- 
ing the atomic frequency and inducing some detuning 
(of the order of lOOg in our simulation), the weakening 
of the blockade effect results in the probability of excit- 
ing to the second manifold to become increasingly strong. 
This means that the polaritonic excitations (or better the 
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photon like polaritonic particles now as we are in the 
dispersive regime) can hop together in bunches of two or 
more from cavity to c&Yity {superfluid regime). 

Simulating spin models: We will now show that in the 
Mott regime the system simulates a XY spin model with 
the presence and absence of polaritons corresponding to 
spin up and down. Let us assume we initially populate 
the lattice only with polaritons of energy ojq — g, in the 
limit (jjd ~ 0^0, Eqs. (2)-(4) becomes 
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(5) 
(6) 



where -P^^''^ — Pi^'^''^ is the polaritonic operator creat- 
ing excitations to the first polaritonic manifold (Fig. 1). 
In the rotating wave approximation, Eq. [7] reads (in the 



H.C. 



interaction picture). Hj — AJ2k=i Pi Pi+i 
In deriving the above, the logic is two step. Firstly 
note that the terms of the type Pi~^^ PiXi^ which inter- 
convert between polaritons, are fast rotating and they 
vanish. Secondly, if we create only the polaritons Pj, ■'^ 
in the lattice initially with energy ujQ — g, then the polari- 
tons corresponding to i^^^^^ will never even be created, 
as the inter-convcrting terms are vanishing. Thus the 



term P, 



(+)t 



can also be omitted. Note that because 
the double occupancy of the sites is prohibited, one can 



(-)t 



identify 
are standard Pauli operators. 



with 



fc' 



where and (j\ 



comes Hi = o'ta 



fe+1 



Then the Hamiltonian be- 
The latter is the standard 



XY model of interacting spins with spin up/down corre- 
sponding to the presence/ absence of a polariton. Note 
that although this is different with optical lattice real- 
izations of spin models, where instead, the internal levels 
of a two level atom are used for the two qubit states 
the measurement could be done using very similar 
atomic state measurement techniques (utilizing the ad- 
vantage of larger distances between sites as well). Some 
simple applications of XY spin chains in quantum infor- 
mation processing such as quantum information transfer 
0, m, [s^l can thus be readily implemented in our sys- 
tem. Very recently a novel idea for efficient cluster state 
quantum computation was proposed in this system where 
the database search and factoring quantum algorithms 
could be implemented using just two rows of cavities psj. 

In conclusion we showed that a range of many body 
system effects, such Mott transitions for polaritonic par- 
ticles obeying mixed statistics could be observed in opti- 
cal systems of individual addresable coupled cavity arrays 



interacting with two level systems. We also proposed pos- 
sible implementations using photonic crystals, toroidal 
microcavities and superconducting systems. Finally we 
discussed the capability and advantages of simulating XY 
spin models using our scheme and noted the ability of 
these arrays to simulate arbitrary quantum networks . 
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